Spin Dynamics in the Magnetic Chain Arrays of SruCi^O^: a Neutron Inelastic 

Scattering Investigation 
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Below T~150 K, the spin arrangement in the chain arrays of Sri4C2404i is shown to develop in two 
dimensions (2D). Both the correlations and the dispersion of the observed elementary excitations 
agree well with a model of interacting antiferromagnetic dimers. Along the chains, the intra- and 
inter-dimer distances are equal to 2 and «3 times the distance (c) between neighboring Cu ions. 
While the intra-dimer coupling is J ~ 10 meV, the inter-dimer couplings along and between the 
chains are of comparable strength, Jm ~ —1.1 meV and J± ~ 1.7 meV, respectively. This remarkable 
2D arrangement satisfies the formal Cu valence of the undoped compound. Our data suggest also 
that it is associated with a relative sliding of one chain with respect to the next one, which, as T 
decreases, develops in the chain direction. A qualitative analysis shows that nearest inter-dimer 
spin correlations are ferromagnetic, which, in such a 2D structure, could well result from frustration 
effects. 



An increasing interest is presently devoted to the com- 
pound Sri4Cu2404i as superconductivity can be ob- 
tained in Ca doped materials under pressure [Q. This 
compound is also a remarkable spin system, which is 
made of two distinct magnetic structures || . One struc- 
ture consists of an array of quantum (s = 1/2) two-leg 
CU2O3 ladders and the other of an array of CuC>2 quan- 
tum spin chains. These two mixed spin subsystems are 
expected to play a crucial role in the occurrence of super- 
conductivity. An important and first question concerns 
the charge transfer process which takes place upon dop- 
ing. Due to the formal Cu valence (+2.5), one expects 
holes to be already present in stoichiometric compounds. 
Recently, it has been proposed that, in such undoped 
materials, most of the holes are located in the chains || . 
This should result in a specific spin distribution in the 
chains. It is the purpose of the present neutron inelastic 
scattering (NIS) investigation to explore the actual spin 
arrangement in the chain arrays of the undoped com- 
pound. A good knowledge of the initial hole distribu- 
tion is definitely required before any investigation of the 
hole-doping process be further developed. A few neutron 
investigations on undoped Sri4Cu2404ihave been previ- 
ously performed. The first measurements were carried 
out on a powder sample [Q. The presence of magnetic 
excitations in the vicinity of 11 meV was confirmed, in 
agreement with susceptibility measurements ||. They 
were attributed to the formation of antiferromagnetic 
(AF) dimers in the chains. The first measurements on 
a single crystal have revealed the presence of several ex- 
citation branches in the same energy range ||. Their 
dispersion, which was determined in a limited part of 



the Brillouin zone, provides evidence of interactions be- 
tween the AF dimers in the chain direction. It was then 
proposed that the dimers are formed between spins sep- 
arated by "2 and 4" times the distance (c) between the 
nearest- neighbor Cu ions in the chains. In these measure- 
ments, however, the chosen scattering plane prevented 
any direct exploration of the chain arrays (fig. la). More 
recently, using a set of eight single crystals, new mea- 
surements were performed H . The dispersion along the 
chains was determined with an improved accuracy. The 
data were shown to compare well to a one-dimensional 
(ID) model of interacting AF dimers (also discussed in 
]To| ) , where the intra and inter-dimer distances (d and 
d'), respectively) correspond to 2 or 3 times the distance 
c. Due to averaging effects inherent to the experimental 
procedure, again no information relative to any trans- 
verse directions could be obtained. In the present inves- 
tigation, a unique single crystal was used and the mea- 
surements were performed on triple-axis spectrometers 
(IN8 and 1T1 at the Institut Laue Langevin in Grenoble 
and at the Laboratoire Leon Brillouin in Saclay, respec- 
tively). The orientation of the crystal was chosen in such 
a way that a rather exhaustive investigation of the chain 
arrays could be developed. 

Our Sri4Cu2404i sample was cut from a crystal sev- 
eral cm long grown by the travelling solvent zone method 
under a pressure of 3 bar oxygen atmosphere [pd| . Its vol- 
ume was about 5 x 5 x 30 mm 3 with the c axis in the 
longer dimension. The crystallographic structure of this 
material can be described with two distinct unitary sub- 
cells, one referring to the ladders, the other to the Cu02 
magnetic chains. For the latter subcell - to be considered 
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hereafter - the space group and the lattice constants are 
at room temperature: Amma, a = 11.456 A, b — 13.361 
A, c = 2.749 A, 0. A sketch of the chain arrays in 
the (a,c) plane is shown in fig. la. In the a direction, 
one observes a doubling of the unit cell (the lattice pa- 
rameter is twice the distance between chains) which is 
due to the relative shift (in the c direction) of one chain 
with respect to the next one Q|J . This relative displace- 
ment is small at room temperature, r\ w 0.17c |Q. Most 
of our neutron measurements were performed with the 
scattering wavevector Q lying in the (a,c) plane. For 
inelastic scatterings, a horizontal collimation of 50'-50'- 
open-open was used, with a final wavevector ky = 2.662 
A . Pyrolytic graphite single crystals were mounted as 
monochromator and analyser, and a pyrolytic graphite 
filter after the sample was used to eliminate higher-order 
flux contamination. Most of the energy scans reported 
here were performed in the range 2 — 20 meV with an 
energy resolution of the order of 1 meV (full width half 
maximum). The temperature could be varied from 1.5 to 
300 K. In general, as shown in fig. 2a, two (and only two 
) well-defined magnetic excitations are observed in the 
(a,c) plane at low temperature. For a few wavevectors, 
however, only a single excitation is detected. It is the case 
for the wavevectors Q = (2.5, 0, Ql), where the unique 
line results from the superposition of two excitations (cf. 
fig. 2a which shows the single peak at Ql = 0.25). More 
surprisingly, it is also the case for the specific wavevectors 
(Qh , 0, 0.65) (see fig. 2b) where an "apparent" extinction 
of one of the two modes occurs. Figs. 3a and 3b display 
the dispersions of the modes observed along the a and 
c directions, respectively. The former figure reports the 
data obtained at Qh — 2 and Qh = 3. For these two 
values of Qh, the data superpose very well yielding an 
excellent definition of the periodicity of the propagation 
along c. The latter figure displays a new and impor- 
tant result. Dispersive modes are also observed in the 
a direction. This result establishes unambiguously that 
a magnetic coupling does exist between the chains. As 
there are 2 non-equivalent neighboring chains in the a 
direction (they form a two-Bravais lattice ) , one is led to 
analyse the two observed branches (E' and E" in fig. 3b) 
as resulting from an unique excitation, which should be 
twice degenerated without such interchain interactions. 
The degeneracy is lifted by the interchain couplings and 
two propagative dispersions are observed. One branch 
is to be considered as the "image" of the other, shifted 
by one reciprocal lattice unit (r.l.u.) 2-n/a (see discus- 
sion below). In contrast to what is proposed in H, the 
two branches do not result from the presence of possible 
magnetic anisotropies. On the contrary, the spin sys- 
tem describing the chains is to be considered as isotropic 
and, accordingly, the two branches are expected to be 
associated with a S — 1 magnetic state. Under the effect 
of a magnetic field, they undergo an identical Zeeman 
splitting giving rise to a unique electron spin resonance 



(ESR) line as observed experimentally With such 
an analysis, no contradiction is found between neutron 
and ESR data. Another striking result revealed by our 
measurements concerns the integrated intensity Iq of the 
observed excitations. The values of Iq measured in the 
chain direction are reported in fig. 4. Several features are 
worth being noted. In particular, extinction is seen to oc- 
cur at specific values of Ql'. 0, 1/2 and 1. For instance, 
for Qh — 3 (fig. 4a), a double peak structure is seen in 
the first half of the Brillouin zone and a single peak struc- 
ture in the second half. All these features must be taken 
into account in any analysis of the magnetic properties 
of the chains. 

As a starting point of our analysis, we consider the 
ID model of interacting AF dimers previously proposed 
|,@. In such a case, we may refer to the description of 
interacting excitons given in When the intra-dimer 
coupling J is much larger than the inter-dimer one J\\ , 
the energy dispersion along the chain is given by 
E/J w (1 - ajj/4) - ct|| (1 + a||)cos(TO27r(3L) + 
(a|/4) cos(to47tQl) « 1 — «|| cos(m27rQi) 
with a; || = J\\/2J. In that expression, Qh is expressed in 
the r.l.u. 2-7r/c, the inter-dimer distance is fixed to d = 2c 
(see discussion below) and the parameter m [— (d+d')/2] 
allows a determination of the inter-dimer distance d! rel- 
atively to d. A fitting procedure where J, a\i and m 
are adjustable parameters is applied successively to the 
two sets of data in fig. 3. Within the experimental 
accuracy, identical values are obtained for ay and to: 
ot\\ = -0.065 ± 0.007 and to = 4.95 ± 0.3. Since a\\ < 0, 
the inter-dimer coupling J|| is deduced to be ferromag- 
netic (F) in that description. The value obtained for 
to suggests the inter-dimer distance d' to be practically 
equal to d' = 3c HJl(J. In that ID description, different 
values are obtained for J : J\ « 12.0 meV and J2 ~ 10.3 
meV for the upper and lower branches, respectively. As 
discussed above, however, the energy difference between 
the two branches results from the small dispersion ob- 
served in the transverse direction. Extrapolating the ID 
description used above to the presence of small interchain 
couplings ( Jj_ <C J), we can write the dispersion observed 
in the a direction as Ej J w 1 — a^cos(7rQ# + (p) with 
a± = J/2 and where the phase factor <p (= or ir) ac- 
counts for the two transverse dispersive branches E' and 
E" . Finally, after fitting the data in fig. 3b (the full and 
dashed lines with | a± |= 0.075 ± 0.007), one obtains the 
following set of values: J « 10 meV, J|| ~ —1.1 meV and 
J± w 1.7 meV. It is established here that the couplings 
between AF dimers along c and along a are of compara- 
ble strength. 

In a second step of our analysis, we consider the data of 
the integrated intensity Iq reported in fig. 4. At low tem- 
perature, as higher energy excitations can be reasonably 
ignored, the quantity Iq provides a good evaluation of the 
"structure factor" Sq — £„(s„s„') exp[iQ(n — n')] where 
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(s n Sn'} defines the static 2-spin correlation functions for 
spins located at sites n and n' . Since J\\, J± <C J, one 
expects Sq to be dominated by the intra-dimer corre- 
lations. Assuming the spins of a dimer to be located 
in a chain at sites n — 1 and n + 1, the different intra- 
dimer correlations are about the same in absolute val- 
ues, i.e. (s n ±is n ±i) w| (s n ±is nT i) |. As the excita- 
tion energy (« 11 mcV) is large compared to the dis- 
persion amplitude (« 1 mcV), the correlation length £ 
of that quantum system is short (£ <C c, a) . Accord- 
ingly, as a first approach, we may limit the sum £„ to 
the inter-dimer correlations between nearest neighbors. 
Then, assuming (in the chain direction) (s n ^is n ±4) w 
(s„±is n±6 ) « -(s nT is n ± 6 ) « -(s„ ± is„± 4 ), we obtain 
from the correlations along the chains the expression 
Iq ~ [1 — cos(47rQi)][l + 2eii cos(to27t(3l)] where the 
first factor accounts for the intra-dimer correlations and 
e\\ is the ratio between the inter- and intra-dimer corre- 
lations (e« = (s n _iS„+ 4 )/(s„_iS n _i}). The first factor 
predicts zero intensities at Ql = 0,1/2,1, exactly as it 
is experimentally observed (see fig. 4). This explains the 
value we have chosen above: d = 2c. In the second factor, 
the coefficient £|| is positive (negative) for AF (F) inter- 
dimer correlations. From the analysis of the dispersion 
curve presented above, e\\ is expected to be negative and 
m f« 4.5. For s\\ = —0.2, one obtains the dashed line in 
fig. 4a. Interestingly, the double peak structure observed 
experimentally is reproduced qualitatively. In the whole 
Brillouin zone, however, the curve remains (almost) sym- 
metric with respect to Ql = 1/2. The occurrence of the 
double peak structure depends directly on the values of 
en: we evaluate 0.15 J$| e\\ |;S 0.25. Within the same sim- 
plified description, the effect of spin correlations between 
dimers located in neighboring chains can be described in 
a similar way and one is led to complete the above ex- 
pression as (for simplicity, the magnetic form factor of 
the Cu 2+ ions is ignored) 

Iq w [1 - cos(4ttQ l )][1 + 2£||Cos(m27rQ L )][l + 
2s j_ cos(ttQ h + m'irQ L )] 

where the last factor takes into account the relative dis- 
placement between neighboring chains: r\ = m'c/2. Sim- 
ilar to the previous case, £j_ is the ratio between trans- 
verse inter- and intra-dimer correlations. Depending on 
whether Qh is even or odd, a phase factor is seen to be 
present or not in the cosine function in the last factor 
of Iq . Considering first the Iq data relative to the exci- 
tation branch E' observed at Qh = 3 (open dots in fig. 
3a), a fitting procedure, with e± and m! as adjustable 
parameters (em = 0.2 fixed), yields the continuous curve 
drawn in fig. 4a, and provides the following evaluations: 
e± = +0.36 ±0.06 and m' = 1 ±0.1. The non symmetric 
behavior observed experimentally is now better repro- 
duced in the whole Brillouin zone. The positive value 
obtained for s± (comparable to £||, as expected) predicts 
the interchain correlations to be AF. Remarkable also is 
the value obtained for m' which shows that the relative 



displacement rj between neighboring chains is of the or- 
der of c/2 at low temperature. In order to test further 
these predictions, the same model with all the above pa- 
rameters fixed (e|| = 0.2, e± = +0.36 and m! = 1) is 
applied to the Iq data relative to the other excitation 
branches (figs. 4b and 4c). According to the descrip- 
tion proposed above, the branch E" observed at Qh = 3 
(black dots in fig. 3) can be viewed as the "image" of 
the excitation branch E' at Qh = 2. The correspond- 
ing structure factor is therefore to be described by the 
same expression as before but without the factor ir in 
the last factor of Iq. One obtains the full curve shown 
in fig. 4b (there is no additional adjustment except for 
an amplitude factor). Again the agreement appears to 
be reasonable in the whole Brillouin zone (for compari- 
son, the curve obtained with the factor ir is shown as the 
dashed line). Considering finally the data obtained at 
Qh = 2, the structure factor relative to the branches E' 
and E" (in fig. 3, the open and full squares, respectively) 
are also reasonably well explained by the same descrip- 
tion: the full and dashed curves shown in fig. 4c arc 
obtained with the factor ir included and not included, re- 
spectively. The only noticeable discrepancy occurs for the 
low-energy data (full squares) in fig. 4c near Ql = 0.65. 
Remember (cf. fig. 2b), it is also the Ql region where 
an unexplained "apparent" extinction of one of the two 
branches occurs. 

The relative displacement 77 can be supposed to vary 
with temperature. Such a change in the lattice structure 
could result in a simultaneous change in the spin arrange- 
ment. This seems to be the case as shown by the data 
reported in figs. 5a and 5b. The former figure represents 
the intensity Aq of two nuclear Bragg peaks, observed at 
Q = (4,0,0.3) and Q = (3,0,0.3). For the chains, these 
peaks correspond to the same Ql = 0.3, but to differ- 
ent Qh (= 3 and 4, respectively). Due to the doubling 
of the unit cell in the a direction, Bragg peaks for even 
values of Ql are expected not to be strongly affected by 
a variation of 77. For odd values of Qh, however, Aq 
should decrease with r] (it is equal to zero for 77 = 0) . 
Figure 5a agrees well with that description and confirms 
the picture suggested above, that neighboring chains are 
slowly sliding along c as T decreases. Simultaneously, 
changes in the excitation spectrum are observed (see fig. 
5b). In that figure, the energy of the single excitation 
shown in fig. 2b, is measured at Qh = 3 and Qh = 2. 
Due to the transverse dispersion, two different values are 
obtained at low temperature. As T increases, however, 
they are seen to merge into the same value w 11.3 meV, 
showing that the interchain coupling becomes ineffective 
above T w 150K. 

The purpose of the analysis sketched above is to draw 
attention on the very peculiar spin distribution appear- 
ing in the chains of Sri4Cu2404i. At low temperature 
(T < 10 K), the interactions between the AF dimers 
develop in two dimensions (2D). For the parameters de- 
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termined above (d = 2c, d! « 3c and r] w c/2), different 
spin arrangements can be proposed. However, as shown 
in fig. lb, there exists one distribution which corresponds 
to a compiete aiternation in the spin dimerisation. It 
is remarkable that a similar situation is found in the 
dimerised phase of the spin-Peierls compound CuGeC>3, 
where a splitting in two excitation branches has also been 
recently observed [ fL3[ . For such a symmetric situation, 
however, if one of the inter-dimer couplings (or both of 
them) are assumed to be ferromagnetic, one is inevitably 
led to a misfit in the spin correlations, i.e. in the spin ar- 
rangement, as illustrated in fig. lb. This picture suggests 
that frustration effects (typically AF couplings between 
next neighbor spins) might also play a role, offering then 
the possibility of alternative descriptions. Anyhow, to 
complete the above analysis, more theoretical and/or nu- 
merical studies are definitely required. The effect of Ca 
doping is expected to modify the charge distribution and 
therefore the spin arrangement in the chain arrays. For 
low doping levels (or non-stoichiometric oxygen contents 
for instance (l4| ) the presence of dimers can still be ex- 
pected, but possibly with a smaller average inter-dimer 
distance d' . This could explain the periodicity of 2 and 
4 reported in For stronger doping levels, however, 
more drastic changes should occur. 

One of us (JPB) would like to acknowledge J.Y. Henry, 
T. Ziman, D. Grempel, F. Mila and G. Bouzerar for nu- 
merous and stimulating discussions. 



FIG. 1. a) Schematic representation of a chain array 
in Sri4Cu2404i, with lattice parameters a and c. r\ is 
the relative displacement between neighboring chains; b) 
spin-dimerisation in the alternation model with an illustra- 
tion of the spin correlations (see text). The full and open 
symbols represent Cu and O atoms, respectively 

FIG. 2. Examples of chain excitation modes at low tem- 
perature (T < 4 K): a) energy scans at [Qh, 0, 0.25) and b) 
at (Q H , 0,0.65). 

FIG. 3. a) dispersion in the chain direction c for excitations 
observed at Qh = 2 and Qh = 3 at low temperature (T < 4 
K); b) dispersion of the two excitations branches E' and E" 
observed in the transverse direction a (T < 4 K). The lines 
are fitted curves (see text). 



FIG. 4. Structure factors of the chain excitations as a func- 
tion of Ql- a) for the branch E' probed at Qh = 3; b) for 
the branch E" probed at Qh = 3; for the branches E' and E" 
(open and full symbols, respectively) probed at Qh = 2. The 
lines are explained in the text. 

FIG. 5. a) Intensity of the two nuclear Bragg peaks ob- 
served at (4, 0, 0.3) and (3, 0, 0.3) as a function of temperature; 
b) energy of the two excitation modes observed at (3, 0, 0.65) 
and (2, 0, 0.65) as a function of temperature. The lines are 
guides to the eye. 
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